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Abstract 

The ordering of scalar fields after a phase transition in which a group G of global 
symmetries is spontaneously broken to a subgroup H provides a possible explanation for 
the origin of structure in the universe, as well as leading to observable effects in condensed 
matter systems. The held dynamics can depend in principle on the geometry and topol¬ 
ogy of the vacuum manifold G/H] for example, texture conhgurations which collapse and 
unwind will exist if the third homotopy group ti^{G/H) is nontrivial. We numerically 
simulate the evolution of texture-like conhgurations in a number of diherent models, in 
order to determine the extent to which the geometry and topology of the vacuum man¬ 
ifold inhuences the held evolution. We hnd that the dynamics is ahected by whether or 
not the theory supports strings or monopoles [characterized by 'Ki{G/H) and 7r2{G/H), 
respectively]. In some of the theories studied, conhgurations with initially spherically sym¬ 
metric energy densities are unstable to nonspherical collapse; these theories are also found 
to nucleate defects during the collapse. Models that do not support monopoles or strings 
behave similarly to each other, regardless of the specihc vacuum manifold. 
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I. Introduction 


Our understanding of large-scale inhomogeneities in the universe, both directly through 
surveys of the distances and velocities of galaxies and indirectly through observations of 
temperature fluctuations in the cosmic microwave background, has advanced signihcantly 
in recent years. This growth in empirical knowledge has necessitated greater precision in 
the extraction of predictions from theoretical models; claiming that a theory is consistent 
with observation requires more care now than it did a decade ago. As a result, it has 
become important to pay close attention to the differences between various manifestations 
of any general scenario. This paper takes up this task for the texture model of structure 
formation, and studies the extent to which the geometry and topology of the vacuum 
manifold in a scalar held theory can exert a cosmologically signihcant inhuence on the 
dynamics of the helds. Although our discussion is phrased in terms of cosmological effects, 
analogous considerations hold for condensed matter systems with spontaneously broken 
global symmetries. 

In the texture scenario, we consider a set of scalar helds $ which transform under a 
global symmetry group G. If the potential V ($) is such that this symmetry is sponta¬ 
neously broken to a subgroup H, the potential will possess a set of degenerate minima. 
This set of held values is known as the vacuum manifold A4, and is isomorphic to the 
quotient space G/H. Typically, at high temperatures the expectation value of the held is 
the same everywhere; as the universe cools the helds will relax to the vacuum manifold, 
and according to the Kibble mechanism the points to which the held evolves will be uncor¬ 
related on distances greater than the Hubble distance at that epoch. The helds will tend to 
order themselves, approaching a constant-held conhguration within a causally connected 
region, but as the Hubble distance grows the helds will always be uncorrelated on cosmo¬ 
logical scales. The resulting gradient energy leads to an approximately scale-free spectrum 
of energy density perturbations, which can serve as seeds for large-scale structure. 
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Depending on the topology of the vacnnm manifold, the scalar helds may be forced to 
leave M in order to smooth themselves ont. This can be seen by considering a conhgnration 
which is conhned to M, and is set to some hxed valne ontside a certain radins. This is 
mathematically eqnivalent to compactifying three-dimensional space to a three-sphere, 
and the resnlting held conhgnration dehnes a map ^ M.. Snch maps are classihed 
topologically by the third homotopy gronp if the conhgnration corresponds to 

a nontrivial element of then it cannot be smoothly deformed to a constant-held 

conhgnration. In this case, it is energetically favorable for the conhgnration to shrink in 
size nntil the energy density reaches that needed to leave the vacnnm manifold. The held 
then “nnwinds” by climbing over the energy barrier to a topologically trivial conhgnration. 
Tnrok [1] snggested that this process of collapse and nnwinding conld lead to seeds for the 
formation of large-scale strnctnre. It was later realized [2,3] that scalar held gradients 
conld lead to density pertnrbations even if were trivial. 

Even thongh density pertnrbations will be prodnced for any theory with spontaneonsly 
broken exact global symmetries, it is reasonable to imagine that the detailed dynamics of 
the helds (and hence the specihc prediction for pertnrbations) will depend on the topol¬ 
ogy of the vacnnm manifold. Indeed, two diherent models with nontrivial which 

snpport textnre conhgnrations with nonzero winding, may nevertheless predict diherent 
evolntions for the collapsing textnre, and snch diherences may manifest themselves cosmo- 
logically. At the same time, it may be possible for conhgnrations which are topologically 
trivial to evolve in ways similar to what we wonld call a textnre. 

In this paper we explore some of these issnes. We consider theories of scalar helds 
transforming nnder global symmetry gronps; the theories are specihed by the gronp G, the 
representation of G nnder which the helds $ transform, and the potential energy E($). Li 
[4] has considered a nnmber of snch theories, and determined what the nnbroken symmetry 
gronp H will be in each model. Bryan, Carroll and Pyne [5] have examined the topology of 
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the resulting vacuum manifolds M = G/H, and calculated the homotopy groups 7ri(7W), 
and (Just as characterizes textures, characterizes cosmic 

strings and Ti 2 {M.) characterizes monopoles.) Here we have chosen eight different mod¬ 
els, listed in Table One, which represent a variety of different held contents and vacuum 
manifold topologies. 

To understand in full detail the cosmological effects of each of these theories, it would 
be necessary to do a number of full-scale numerical simulations of large-scale structure 
formation or CMB anisotropy in each model, such as those which have been performed for 
the models with M. = [6,3]. We have adopted a more modest approach, since our goal 

is simply to determine whether the choice of vacuum manifold exerts an inhuence over 
the dynamics. Our strategy is therefore to set up comparable held conhgurations in the 
various theories, each corresponding to a single texture. We then follow numerically the 
evolution of the conhgurations as they collapse and evaporate away, keeping track of the 
distribution of energy into potential, kinetic and gradient energy, as well as the asymmetry 
of the collapse as measured by the quadrupole moments of the total energy density. (For 
analogous studies of the usual texture models, see [7,8].) In each model we performed three 
simulations: one in which the energy density was initially spherically symmetric, one which 
was initially prolate (cigar-shaped), and one which was initially oblate (pancake-shaped). 

As detailed below, we hnd that the models studied fall naturally into two diherent 
categories of behavior. Within each of the two classes, the evolution of the energy densities 
and quadrupole moments is relatively similar for each type of collapse, while there is a 
noticeable difference between the two categories. In the eight models we studied, the 
division corresponds precisely to whether the theory supports strings or monopoles; be., 
to whether or not and Ti 2 {M.) are trivial. All of the models without these defects 

evolved in one way, while those which do support defects evolved in a distinct fashion. 
Upon examination of the simulations, we discover that in some of the models defects are 
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nucleated in the process of collapse; however this is not always true, and is therefore not 
a sufficient explanation for the division into two types of collapse. 

II. Initial Configurations 

The models we consider are described by Lagrangians of the form 

T = -]/($) . (2.1) 

Here, $ stands for a collection of N scalar helds which transform under some representation 
of a global symmetry group G. The product ■ $2 is the invariant inner product appro¬ 
priate to the representation; for a vector we have $1 ■ $2 = Xla $i *$2 (where an asterix 
denotes complex conjugation), while for a matrix we have $1 ■ $2 = ~ 

Tr($|$ 2 )- The potential H($) consists of terms quadratic and quartic in the helds. 

We are interested in constructing, for each theory, initial held conhgurations which 
are wholly in the vacuum manifold and which have unit winding number. We follow the 
procedure outlined in [5]. The basis of this procedure is the fact that every held value in 
the vacuum manifold may be obtained by starting with a hxed vacuum expectation value 
($) and acting on it by an element of the symmetry group G. The initial conhguration 
can therefore be written 

l>(x)=Mx)W. (2.2) 

where /u(x) is a position-dependent element of G. (The notation //(x)($) indicates the 
action of n on the representation carried by ($); for vector representations it will be matrix 
multiplication, while for matrix representations it will be conjugation.) By considering 
transformations /u(x) which go to the identity at spatial inhnity, n dehnes a map n : ^ 

G, and hence an element of 7r3(G). Elements of ti^{G/H) are related to those of 7r3(G) by 
the exact homotopy sequence 

7r3(G) ^ ti^{G/H) ^ Ti2{H) . (2.3) 
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Since ti 2 {H) = 0 for any H, exactness implies that every element of ti^{G/H) is the image 
nnder (3 of some element of 7r3(G). To constrnct a conhgnration with winding nnmber 
one in the vacnnm manifold, it is therefore necessary to hnd a map // with an appropriate 
winding nnmber in 'n's{G). 

The map fi can always be written as a composition, 


H = a o fj, , 


(2.4) 


where /i is a hxed map from to SU(2) and a maps SU(2) to G. If we choose Jl to 
have winding nnmber one, then fi will represent the same winding nnmber as a. For 
Jl : ^ SU(2) we choose 


^(x) 


f c^+ izs^ {ix + y)^x 
\{ix-y)s^ c^-izs^ 


(2.5) 


Here we have introdnced the notation = cos(x), — sin(x), x(’') a fnnction of r = 
with x(0) = 0 and x(oo) = ^md x = x/r. Constrnction of an appropriate 
held conhgnration then comes down to choosing the map a : SU(2) — G. 


We begin with a theory in which S0(4) breaks down to S0(3). In this theory the 
vacnnm manifold is SO(4)/SO(3) = S^; it is the simplest theory containing textnres, and 
may be thonght of as the standard against which other models shonld be compared. The 
held content consists of fonr real scalars (a = 1,...4), transforming nnder a vector 
representation of S0(4). The action of an S0(4) matrix O^b is thns h-T he 
potential is given by 

H($“) = -^2)2 . (2.6) 


This potential is minimized when the helds attain a vev of the form 


($) = 


^ 0 ^ 

0 

Vo/ 


(2.7) 
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breaking SO(4) to an SO(3) snbgronp, namely the snbgronp consisting of 4 x 4 matrices 
with a one in the npper left corner, zeroes in the rest of the hrst row and colnmn, and SO (3) 
matrices in the lower right 3x3 snbblock. The map a mnst now be a homomorphism from 
SU(2) to SO(4) with winding nnmber one. [In fact 7r3(SO(4)) = Z © Z, so the concept of 
winding nnmber is not nniqnely dehned. The map we choose is that which indnces the 
generator of 7r3(SO(4)/SO(3)).] If we write a general element of SU(2) in terms of fonr 
real parameters with ScaCa == 1, 


f 60 + ies 62 + i6i \ 
y -62 + i6i 60 - *63 y ’ 


then an appropriate map is given by 


a(^) 


/ eo 

63 

62 

61 

-63 

60 

-61 

62 

-62 

61 

60 

-63 

V-6i 

-62 

63 

60 


We may think of a as replacing complex entries by real 2x2 matrices: 


( 2 . 8 ) 


(2.9) 


a + ib 


a b 
—b a 


The composition of a with /u(x) resnlts in a map /*(x) given by 

/ 


^(x) = 

which leads in tnrn to a held conhgnration 

$(x) = V 




y^x 

xs^ \ 

-ZSx 

^x 

—XSx 

y^x 

-y/x 

XSx 

^x 


—XSx 

~y^x 

ZSx 



( 


\ 


-zs^ 

-yfx 

—xs^ 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


The model in which SO(5) breaks down to SO(4) is closely related to the previons 
one. The helds are arranged into a hve-vector with potential as in (2.6). The vacnnm 
manifold is SO(5)/SO(4) = for which 7r3(S'^) = 0; there are therefore no conhgnrations 
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with nonzero winding nnmber. Nevertheless, we can set np conhgnrations which look jnst 
like those in the previons theory and see how they evolve. This serves as a test of the 
importance of the topology of the vacnnm manifold in textnre scenarios. 


The conhgnrations we consider are maps from to which lie at a constant distance 
from the north pole in 5”^. They are given by 


(2.13) 


$(x) == V —ays^ 




The parameter a determines the distance from the north pole; a = 1 corresponds to the 
eqnator of while a = 0 pnts the entire held conhgnration at the pole (a single point). 

Another closely related theory featnres SO(5) breaking down to SO(3). This model 
provides an interesting comparison with SO(4)/SO(3), since in this case ti^{G/H) is the 
hnite gronp Z 2 . There are ten real scalar helds, arranged into two hve-vectors and $25 
with potential 


V (<!>?, Kl) = Ai (<I>;<I>; -v}f + M (<I>‘<I>‘ - vlf + r, (<i>;<i>5)- 


(2.14) 


The helds acqnire vevs of the form 


/ni\ /U 

0 V2 

($1) = 0 , ($2) = 0 

0 0 

Vo/ Vo, 


(2.15) 


We act on these helds by a 5 x 5 matrix with (2.11) as the npper left 4x4 snbblock, a 
one in the lower right corner, and zeroes elsewhere. (Note that a single SO(5) matrix acts 
simnltaneonsly on the two vectors.) The resnlting conhgnration is 


$1 = Vi 


, $ 2—^2 


(2,16) 



We next consider a theory where SU(3) is completely broken {H = 0). (The ap¬ 
pearance of textnres in a particle physics model with this symmetry-breaking pattern was 
examined in [9].) The lower homotopy gronps of the resnlting vacnnm manifold are the 
same as those of SO(4)/SO(3), althongh the nnmber of helds and the dimensionality of 
the vacnnm manifold are both greater. The helds consist of two complex three-vectors 
and $2 5 ^ total of twelve real degrees of freedom. The potential is given by 

2 2 2 

y ($?,$^) = Ai -nf) -f-A2 ($^$2 - + • (2.17) 


The helds attain vacnnm expectation valnes of 


(l-l) = jj j 



( 2 . 18 ) 


In this case the map a : SU(2) —SU(3) is simply inclnsion of a 2 x 2 matrix into the 
npper left corner of a 3 x 3 matrix. Performing the transformation ^ = a oji on the two 


vectors then yields 

/ c^ + izs^ \ 
$1 = {ix - y)s^ 

\ 0 / 


$2 


/ (tx+ 

V2 I -^ZS^ 1 


(2.19) 


The theory in which SO(3) breaks down to SO(2) is once again similar to SO(4) 
breaking to SO(3). The resnlting textnres are known as “Hopf textnres” [10]. The helds 
comprise a real three-vector $“ with potential of the form (2.6), which leads to a vev 


($) = 



( 2 . 20 ) 


The vacnnm manifold is SO(3)/SO(2) = for which both 7r2 and tts are Z. The theory 

therefore contains monopoles as well as textnres. Unlike the two cases above the map a of 

(2.9) will clearly not snfhce. Instead with g as in (2.8), we set 

! el + el - el-el 2(6162 + 6063) 2(6163-6062) \ 

Q;( 5 r) = 2(6162-6063) el - el el-el 2(6263 + 6961) . ( 2 . 21 ) 

\ 2(6163 + 6062) 2(6263-6061) el - el-el + el ) 
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This map is the standard double cover of SO(3) by SU(2). (The winding number is 
nevertheless one; there are no single covers.) Its composition with /u(x) from (2.5) results 


in 

( c^ + { 2 x‘^ — l)s^ 2 {xys'^ — zc^s^) 
2 {xys‘l + zc^s^) c\ + (2y2 _ 1)^2 
2{xzsl,-yc^s^) 2{yzs'^^ +xc^s^) 

The resulting held conhguration is straightforward: 


2{xzs'^ +yc^s^) 
2{yzsl 




+ (2^2 - l)s^ 


/4 + ( 2 x 2 _i) 4 \ 
$(x) = V 2{xys\ + zc^s^) 
\2{xzsl-yc^s^) j 


( 2 . 22 ) 


(2.23) 


Adding an additional three-vector to the above theory yields a model in which SO(3) 
is completely broken. The vacuum manifold is therefore SO(3) itself, which supports Z 2 
strings as well as textures. The potential is of the form (2.14), leading to vevs of the form 



Acting on these helds with (2.22) gives the following held conhguration: 

/ 4 + (2£2 _ i)s2 \ / 2{xys\ - zc^s^) \ 

$i(x) = vi ( 2 {xys\ + zc^s^) , $2(x) =^214 + { 2 f - 1)4 

\2{xzsl-yc^s^) j \2{yzsl + xc^s^) j 


(2.24) 


(2.25) 


We turn now to a model in which an SO(4) symmetry is broken to U(2). This is 
accomplished by a set of six scalars which transform in the antisymmetric tensor repre¬ 
sentation of SO(4), The action of an SO(4) transformation O on the helds is 

given by $ H-> The quartic potential invariant under this symmetry is given by 

t/($“^) = -lm2Tr($^$) + ^Ai[Tr($2)]2 + lA2Tr($^) , (2.26) 

^ hi hi 


where refers to matrix multiplication. The helds acquire a vev of the form 


($) = V 


(° 

1 

0 


-1 

0 

0 

0 

0 

0 

0 

1 

\ 0 

0 

-1 

0 / 


(2.27) 
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In this case we cannot use (2.11) to act on the vev, since this transformation is actually 
in the unbroken subgroup H = U(2), and hence leaves (2.27) invariant. (The generator of 
7 r 3 (SO( 4 )/U( 2 )) and the generator of 7 r 3 (SO( 4 )/SO( 3 )) are inherited from the two different 
generators of 7 r 3 (SO( 4 )) = Z © Z.) Instead we take the SO(3) transformation (2.22) and 


include it in SO (4) by putting a one in the lower right corner and zeroes elsewhere. Acting 
the result on (2.27) yields a held conhguration 



0 

4 + (2h2 - 1)4 

-2{xc^s^ + yzs"^^) 

-c2 - 

(2h2 - 1)4 

0 

2(xzs^ -yc^s^) 

2{xc^ 

sx + y^^l) 

-2{xzsl. - yc^s^) 

0 

-2{xzs\ - yc^s^) 

-2{yzs‘^^ + xc^s^) 

-4 - (2^2 -1)4 


2{xzs‘l - yc^s^) \ 
2(yzs| + xc^s^) 

4 + ( 2^2 - 1)4 
0 / 
(2.28) 


Our hnal model features SO(5) breaking down to SO(3) x SO(2) x Z 2 . There are 
fourteen scalar helds transforming in the symmetric and traceless tensor representation of 
SO(5): Tr$ = 0. The potential is 


Tr($2) + lAi[Tr($2)]2 + 4^ Tr($4) , 

^ hi hi 

leading to a vev of the form 



0 

0 

0 


0 

2 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

-3 

0 

\0 

0 

0 

0 

-3/ 


(2.29) 


(2.30) 


The map /u(x) we have already found in (2.11). Its action on the vev yields the held 
conhguration $ = which in components is 


$ = n 


/2 — 5x^s'^ 

-5xys^ 

bxzs'^ 

5xc^s^ 


-5xys^ 

2 - bfsl 

5yzs‘l 

2 - 5z^sl 

^y^x^x 

0 

5xzs^ 

5yzs‘^ 

—5zc^s^ 

0 

5xc^s^ 


—5zc^s^ 

2-54 

0 

\ 0 

0 

0 

0 

-3/ 


(2.31) 


To compare the diherent theories, we would like to compare the evolution of conhgura- 
tions which are in some sense “the same,” despite being in diherent models. For example. 
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it would be reasonable to compare the collapse of initially spherically symmetric textures. 
In all of the above held conhgurations, the initial energy density is spherically symmetric 
if we choose the radial function 


x(r) = 2 arctan(r) . 


(2.32) 


Indeed, with this choice the energy density prohles are identical (up to an overall normal¬ 
ization); they are given by 


p{x,y,z) oc 


1 

(l + r2)2 • 


(2.33) 


In this sense, the initial conhgurations we consider are directly comparable. 


However, there is a more restrictive dehnition of spherical symmetry: that a conhg- 
uration which is rotated about the origin in space can always be brought to its original 
form by a global symmetry transformation. Not all of the theories we consider allow 
conhgurations which are symmetric in this strict sense. It is straightforward to check, 
by considering inhnitesimal spatial rotations and transformations in the symmetry group 
G, that all of the textures in the models without defects [SO(4)/SO(3), SO(5)/SO(4), 
SO(5)/SO(3), and SU(3)/0] are spherically symmetric, and in addition the conhguration 
inSO(5)/[SO(3) X SO(2) X Z 2 ] is as well. In the remaining theories it is impossible to set up 
a truly spherically symmetric conhguration. For the theories based on SO(3)/SO(2) = S'^ 
and SO(4)/U(2) = Rp2^ the vacuum manifolds are two-dimensional (and the map from 
space to the vacuum manifold is onto); therefore, the preimage of any given point in 
the vacuum manifold will be a one-dimensional region of space, which in general will not 
be taken into another such one-dimensional preimage under a spatial rotation. In the 
SO(3)/0 = theory, meanwhile, the spatial rotation has precisely two hxed points (the 
origin and the point at inhnity), while the group action has none; therefore, a rotation 
cannot be undone by a group transformation. 

As we shall see, this lack of true spherical symmetry will manifest itself in our simula¬ 
tions. This can be traced to the following effect: if a conhguration is not truly symmetric 
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but has a symmetric energy density, deformations of the conhguration along different coor¬ 
dinate axes will not yield conhgurations of equal energy. Hence, there will be an instability 
to non-spherical collapse. We will see this explicitly in our simulations of textures with 
initially spherically symmetric energy densities. Furthermore, we also perform simulations 
of oblate and prolate textures, obtained from those with symmetric energy densities by 
deforming along one axis. For the truly symmetric initial conhgurations, the axis chosen is 
of no consequence, and the energy density of the deformed texture is simply the deformed 
energy density of the symmetric texture. The conhgurations which are not truly symmet¬ 
ric, however, have a preferred direction [the x-axis for (2.23), and the z-axis for (2.25) and 
(2.28), with our conventions]. A deformation orthogonal to the preferred direction results 
in a triaxial energy density with no rotational symmetry, while a deformation along the 
preferred direction yields an axially symmetric energy density. In SO(3)/0 this axially sym¬ 
metric pattern is very similar to that of the models which allow spherical symmetry, while 
in SO(3)/SO(2) and SO(4)/U(2) the resulting energy density is toroidal. In Figure One 
we have plotted the energy densities in two-dimensional slices through two representative 
deformed conhgurations, SO(4)/SO(3) and SO(3)/SO(2). 

The energy density prohles of the oblate and prolate conhgurations are therefore 
not identical in diherent theories, and consequently the comparison between the models 
is not as direct as in the initially symmetric case. We will see below, however, that 
the diherence is not crucial. For our simulations of deformed textures, we have always 
performed the deformation along the preferred axis, so that the comparison is always 
between conhgurations with axially symmetric energy densities. 

III. Numerical Implementation 

We integrate the held equations of motion using a standard staggered leapfrog algo¬ 
rithm, using a method that is second-order in time but which uses fourth-order spatial 
diherences (as in [8], except we have added fourth-order spatial diherences for better res- 
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olution). The field equations for the above Lagrangian with a general potential 1^($) 


are 

-+ —=0. (3.1) 

dt^ ^ ’ 

Here, the field $ is vector or tensor valued, depending on which theory we are integrating. 
In the staggered leapfrog method, the momenta of the field tt = $ are defined on half¬ 
integer timesteps. Using finite differences as an approximation to the differential operators, 
the expression for the field at discrete timestep n -|-1 and Cartesian grid location x = iAx, 
y = jAx and ^ = kAx is 

= A* (3.2) 


and from the field equations, combined with the momentum definition, we obtain the 
expression 


n+l/2 _ ri-1/2 

‘^ijk ~ '^ijk 


+ ^lj+2,k + ^"j,k+2 - 16$r+l,,fc - 16$-, + 




+ ^1^-2,jk + ^"1,3-2,k + ^'ij,k-2) 


At 


dV 


(3.3) 


for TT at timestep n + 1/2. 


We use Neumann boundary conditions. In our simulations, once the initial texture 
configuration collapses, massless and massive radiation hits the wall, then some of it is re¬ 
flected and influences the subsequent evolution of the field inside the box. In our production 
simulations, we take the radius of the texture to be 10 gridzones, while the simulation vol¬ 
ume is 64 X 64 X 64 gridzones. With this geometry, reflected radiation does not influence 
the texture configuration until after it has collapsed and the decay products (either defects 
or radiation) have had time to propagate away from the texture core to a distance of order 
the initial size. 


The initial texture configurations are discussed above for the various theories which 
we have simulated. For the purposes of the numerical simulation it is convenient to have an 
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initial condition with zero energy at the bonndaries of the simnlation volnme; we therefore 
modify the radial fnnction x(r) given in (3.32) so that the helds are initially constant 
ontside a hxed radins: 


X(r) 


2 [l — arctan(r) for r < i? , 

TT for r > i? , 


(3.4) 


with (3 and R set snch that x goes smoothly to tt at the hnite radins R. We also examined 
textnres in a theory where SU(3) was spontaneonsly broken to SO(3); in this model, 
however, no natnral choice of x(?’) led to a symmetric energy density. We have therefore 
not inclnded this model in onr comparisons. 

There are two conditions on the timestep for a global symmetry breaking theory. 
One is the nsnal Conrant condition \/2At < Ax. In addition, the sonrce term from 
the derivative of the potential dV/d^ adds an effective mass term to the eqnations. For 
stability, the timestep mnst be short enongh that variations in the held are well resolved. 
One way to see this is to look at a vector theory in the sitnation where the held is spatially 
constant. We then have an eqnation 


$ = - A (|$|2 ^ 


(3.5) 


where the ehective mass is A(|$p — n^). If we normalize v to be 1, then when $ is close to 
zero and at the core of the potential, the ehective mass is at a maximnm, and the solntion 
to the above diherential eqnation has a wavennmber of To resolve the solntion for 

this wavennmber, we typically take the timestep to be 1/10 of the wavelength. 

For the parameters of the scalar potential, we take the qnartic self-conpling of the 
held to be A = 0.1, and the magnitnde of the vacnnm expectation valne to be n = 0.1. 
The SO(5)/SO(4) model contains the additional parameter a, which characterizes where 
on the vacnnm manifold the initial conhgnration lies; we have nsed a = 0.1 in the 
simnlations shown. (Additional simnlations with other valnes of a led to indistingnishable 
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results.) Together with a timestep At ~ 0.1 our code is stable for all collapses in the 
theories studied. 

To ensure stability and the accuracy of the code, we check that the total energy within 
the simulation volume is constant. We hnd that for the above parameters, the energy is 
typically constant to within 1%. There is often a slight increase due to the production 
of massive radiation. This is due to an increase in the effective mass of the held, and 
subsequently a poorer resolution of huctuations in the held (discussed above) in time. 

IV. Results 

For each of the eight models listed in Table One, we performed three simulations of 
single-texture collapse; one with a spherically symmetric initial conhguration, one with an 
oblate initial conhguration, and one with a prolate initial conhguration. The oblate and 
prolate conhgurations were obtained by scaling the conhguration along either one or two 
axes by a factor of 0.5. At each timestep in the simulations we calculate the total energy E, 
as well as the separate values of the kinetic energy gradient energy Eg^ and potential 
energy Ey- 

E =Ek + Eg -|- Ey 



To measure the asymmetries in the collapse, we computed the quadrupole moments of the 
total energy in a sphere of radius the initial texture size: 

Qij = f d^x {3x'^x^ — (4.2) 

J r<R 

where p is the total energy density. Since we set up our conhgurations with initial asym¬ 
metries along the coordinate axes, we are interested in the moments Qu which measure 
the asymmetry along the x^ axis. 
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We examine first the four models listed in Table One which do not support strings 
or monopoles; = Ti 2 {M.) = 0. In Figure Two we have plotted the integrated 

energies (total, gradient, kinetic and potential) for the evolution of initially symmetric 
conhgurations in each of these theories. The total energy is normalized to 100 for purposes 
of comparison. These graphs reveal several features which are common to all of the collapse 
simulations we studied. The conhguration starts out with zero kinetic and potential energy. 
As the texture collapses the potential energy grows slightly, but remains a small part of the 
total energy throughout the simulation. The kinetic and gradient energies move toward 
equipartition, although this is not necessarily completely successful. In every simulation 
we performed, the graphs of the kinetic and gradient energies approach each other and 
cross at least once; the point where they are hrst equal serves as a convenient measure of 
the timescale for the collapse. In the simulations shown in Figure Two, equality is reached 
after approximately 75 timesteps. 

A remarkable feature of these four graphs is their similarity with respect to each 
other. In the SO(4)/SO(3) and SO(5)/SO(4) theories the kinetic and gradient energies do 
not equipartition quite as effectively as in the SO(5)/SO(3) and SU(3)/0 theories, but the 
difference is small. As we shall see below, the differences within this set of models are much 
less than the differences between them and the models which support defects (strings or 
monopoles). 

The results of the SO(5)/SO(4) simulation are especially noteworthy, since in that the¬ 
ory the initial conhguration is topologically trivial [7r3(SO(5)/SO(4)) = 0]. Nevertheless, 
the evolution of the energy densities proceeds in a manner practically indistinguishable 
from that in SO(4)/SO(3). (The fact that such conhgurations would collapse was noted 
in [2]; our simulations provide dynamical evidence that the collapse is almost identical to 
those in theories with nontrivial The only noticeable diherence is in the potential 

energy, which remains essentially zero throughout the SO(5)/SO(4) collapse, while rising 
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slightly in the other models. This is due to the fact that, in the theories with topologically 
nontrivial conhgurations, the helds must climb out of the vacuum manifold in order to 
unwind, which necessarily induces a nonzero potential energy. Since the potential energy 
is nevertheless only a small contributor to the total energy of the texture, these simula¬ 
tions indicate that the behavior characteristic of “texture collapse” may arise even without 
nontrivial topology. (Of course the frequency with which such behavior actually does arise 
in a cosmological context is not addressed by our simulations, and may be different in the 
different models.) 

We have not shown the evolution of the quadrupole moments for these simulations, 
since they vanish throughout the course of the collapse. This simply indicates that these 
models remain spherically symmetric if they begin in a symmetric conhguration. 

In Figure Three we have plotted the same quantities as in Figure Two, in the same 
set of models, this time for initially oblate conhgurations. (The results for prolate collapse 
will not be shown, as they are essentially identical to the oblate case.) Once again, the 
results of these four simulations are very similar to one another. They are also similar to 
the results for initially symmetric conhgurations; again we see the equipartition of kinetic 
and gradient energies, while the potential energy is much smaller. The timescale for these 
collapses is slightly less than that in the symmetric case; the kinetic and gradient energies 
hrst become equal after approximately 60 timesteps. This difference can be traced to the 
fact that the initial conhgurations are somewhat smaller, as they are obtained from the 
symmetric conhgurations by shrinking by 50% along the x-axis. 

The quadrupole moments for the evolution of the initially oblate conhgurations in 
models without strings or monopoles are shown in Figure Four. The three quantities Qxx: 
Qyy and Qzz are plotted; only two curves are visible because the initial deformation is along 
the x-axis, and hence Qyy = Qzz throughout the simulation. The quadrupoles begin at a 
certain value, head toward zero, and cross the axis before settling smoothly toward zero. 
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The physical interpretation is that the initially oblate conhgnration moves toward spherical 
symmetry, then overshoots to become somewhat prolate. The qnadrnpole moments are 
only calculated in a sphere of radius equal to the initial texture size, so, as the radiation 
leaves the sphere, the qnadrnpole goes to zero. The remarkable resemblance between the 
different simulations is the strongest evidence that texture collapse is essentially the same 
in these four models; not only do the total energies partition in similar ways, but the 
spatial distribution of the total energy density for nonspherical initial conhgurations also 
evolves similarly. 

We turn next to the four models listed in Table One which support strings or monopoles. 
(The initial conhgurations that we set up are pure texture, since the helds are everywhere 
in the vacuum manifold, but these models allow topologically stable defect conhgurations.) 
In Figure Five we have plotted the integrated energies (total, gradient, kinetic and poten¬ 
tial) for the collapse of conhgurations with initially spherically symmetric energy densities 
in these models. The four graphs are similar to each other, although we shall see below 
that these models are distinguishable by other measures. As in the models that do not 
support defects, the kinetic and gradient energies evolve toward equipartition, while the 
potential energy is always signihcantly smaller. A comparison between Figures Two and 
Five, however, reveals that there are important diherences between the spherical collapses 
in models with and without strings or monopoles. The most signihcant quantihable dif¬ 
ference is in the timescale for collapse; in the models which allow defects, the kinetic and 
gradient energies become equal after only 45 timesteps, in comparison with 75 timesteps 
for the models of Figure Two. 

Another dramatic difference is manifested in Figure Six, which shows the qnadrnpole 
moments for the collapses of initially spherical conhgurations in models which support 
defects. The notable feature of these plots is that the quadrupoles in three of the four 
models [SO(3)/SO(2), SO(3)/0, and SO(4)/U(2)] do not vanish throughout the simula- 


19 



tion, even though the initial energy densities are spherically symmetric. As discussed at 
the end of Section Two, the asymmetries indicated by the nonzero quadrupoles can develop 
because the initial held conhguration is not truly spherically symmetric, even though the 
energy density is. The collapse of these textures is therefore different in an important way 
from the collapse of truly symmetric textures, since the energy densities evolve asymmet¬ 
rically. On the other hand, it is also important to recognize that in the remaining model, 
SO(5)/[SO(3) X SO(2) X Z 2 ], which is truly symmetric and does not develop a nonzero 
quadrupole, the partition of energy densities nevertheless resembles the other models which 
support defects more than those which do not. This is why it is more appropriate to clas¬ 
sify the theories by whether or not they support strings or monopoles, rather than whether 
or not they admit spherically symmetric texture conhgurations. 

Closer examination of the held values during these simulations reveals that the asym¬ 
metry leads to dramatic effects, in the form of nucleation of topological defects. Figure 
Seven is a three-dimensional rendering of contours of total energy density and potential 
energy at a single timestep during the SO(3)/SO(2) simulation. This theory supports 
monopoles as well as textures, and the disk-shaped regions of high potential energy are 
monopoles and antimonopoles. A sequence of such pictures throughout the simulation 
reveals that the texture initially collapses in a nearly symmetric fashion. It then begins to 
deviate from spherical symmetry (while remaining axially symmetric), and a monopole- 
antimonopole pair nucleates at the center. These travel away from each other toward the 
edges of the simulation region, while eventually another pair nucleates at the center. The 
antimonopole of the second pair follows the monopole of the hrst, and vice-versa. Presum¬ 
ably these will annihilate, although the hrst pair leaves the box before this can occur. The 
simulation of SO(4)/U(2) proceeds along the same lines, with identical results. This is not 
surprising, since the vacuum manifold SO(4)/U(2) = is related to the vacuum man¬ 
ifold SO(3)/SO(2) = by identihcation of antipodal points; the local geometry of both 
vacuum manifolds is therefore the same. Laboratory experiments with nematic liquid crys- 
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tals [11] have seen evidence of similar behavior; the vacnnm manifold in that case was RP^, 
and it was fonnd that textnres tended to decay into monopole-antimonopole pairs. Mean¬ 
while, the distinct geometry of SO(3)/0 = leads to dramatically different evolntion; 
in this model there are no monopoles, bnt the textnre collapse leads to formation of a loop 
of cosmic string, as shown in Fignre Eight. This loop grows and intersects the edges of the 
simnlation region, so we cannot follow its entire evolntion. The SO(5)/[SO(3) x SO(2) x Z 2 ] 
theory, althongh admitting both strings and monopoles, collapses symmetrically, and no 
defects are formed. 

Of conrse, in a cosmological context it wonld be very nnlikely for a collapsing textnre 
to begin in a perfectly symmetric conhgnration. To this extent, the instability of some 
theories is of less importance than the evolntion of initial conditions which are noticeably 
asymmetric. Fignre Nine presents the behavior of the integrated energies in the fonr 
models which snpport defects, beginning from oblate initial conditions. Again, the fonr 
plots appear similar to each other, while differing from the analogons plots in Fignre Three 
for models withont strings or monopoles. The timescale for collapse is now approximately 
35 timesteps, and the eqnipartition of kinetic and gradient energies occnrs rapidly and 
effectively. 

Fignre Ten, which plots the qnadrnpole moments for these collapses from oblate initial 
conditions, indicates that the models with defects can differ from each other as well as from 
the models withont defects. (Althongh it is not evident from these plots, we have checked 
that defects are nncleated according to the same pattern which appeared in symmetric 
collapses; monopole/antimonopole pairs are prodnced in SO(3)/SO(2) and SO(4)/U(2), 
a string loop is prodnced in SO(3)/0, and no defects are prodnced in SO(5)/[SO(3) x 
SO(2) X Z 2 ].) The general pattern familiar from Fignre Fonr, where the oblate conhgnration 
became prolate before symmetrizing, is seen here as well. The qnadrnpole moments in the 
SO(3)/SO(2) and SO(4)/U(2) models, whose vacnnm manifolds are locally eqnivalent. 
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are indistinguishable, but these two appear different from SO(3)/0 and SO(5)/[SO(3) x 
SO(2) X Z 2 ], which also appear slightly different from each other. As these differences 
represent distinct evolutions of the spatial distribution of the total energy density, they 
indicate the possibility that the different models could lead to distinguishable cosmological 
density fluctuations. 

V. Discussion 

We have examined the dynamics of texture collapse in a variety of scalar field theories 
with spontaneously broken global symmetries. Our goal was to determine whether the 
geometry and topology of the vacuum manifold would affect the evolution of analogous field 
configurations in different models. The answer is clearly “yes”; however, the cosmological 
significance of the effects observed is still unclear. 

We have been led to two basic results. First, the collapse of a single texture configu¬ 
ration depends on whether the theory under consideration admits the existence of cosmic 
strings or monopoles. Theories which do allow such defects have textures which collapse 
somewhat more rapidly, and in the process of collapse may nucleate loops of string or 
monopole/antimonopole pairs. Second, some theories with nontrivial TTsi-M.) admit tex¬ 
tures with spherically symmetric energy densities, but not texture configurations which 
are strictly spherically symmetric. Such configurations collapse nonspherically. These two 
results are related; it is the configurations without true spherical symmetry that nucleate 
defects during collapse. In the wide class of models considered in [5], all of the theories 
without monopoles or strings either allow spherically symmetric textures or do not have 
textures at all. 

From the point of view of cosmology, texture collapse is not the only important process 
in a theory which admits strings or monopoles; in such models the Kibble mechanism 
predicts that there will be an appreciable number of defects per Hubble volume, and the 
gravitational perturbations produced by these defects will presumably be comparable in 
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importance to those produced by collapsing textures. Indeed, the results of our simulations 
may indicate that textures are even less important in such models than might have been 
naively expected. The laboratory experiments of [11], involving a condensed-matter system 
which admitted strings, monopoles and textures, found that textures were rarely produced 
in a phase transition, and when a texture conhguration was introduced by hand it decayed 
into monopole/antimonopole pairs. Our results, in which textures with initially perfectly 
symmetric energy densities were seen to collapse into defects, may indicate that textures 
will rarely form in certain theories, with defect formation being more likely. Contrariwise, 
the collapse of textures may provide an additional mechanism for defect production. The 
importance of these effects to cosmological structure formation can only be answered by 
large-scale numerical simulations of the evolution of the helds as the universe expands. 

For the models which are free of strings and monopoles, our results provide evidence 
that simulations of any specihc simple model [such as SO(4)/SO(3)] do not miss any 
important physical effects. We hnd that the collapse of textures proceeds in essentially the 
same way in all such theories; to this extent, simulations in one model are relevant to other 
models as well. Of course, there are a number of important issues which our simulations 
do not address, such as the frequency and distribution of texture collapse, which could in 
principle be different in distinct models. 

Finally, our examination of the SO(5)/SO(4) model provides evidence that events re¬ 
sembling texture collapse can occur in the absence of nontrivial topology, as suggested in 
[2]. Indeed, the quantitative characteristics of the collapse in this model are indistinguish¬ 
able from those in the conventional SO(4)/SO(3) model. Again, the relative frequency 
of such events in the different models is not addressed by our simulations. More detailed 
studies will be necessary to understand the precise differences between cosmological models 
based on the theories examined here. 
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G 

H 

N 

7ri(G/H) 

7T2iG/H) 

TisiG/H) 

dim{G/H) 

SO(4) 

SO(3) 

4 

0 

0 

z 

3 

SO(5) 

SO(4) 

5 

0 

0 

0 

4 

SO(5) 

SO(3) 

10 

0 

0 

Z 2 

7 

SU(3) 

0 

12 

0 

0 

z 

8 

SO(3) 

SO(2) 

3 

0 

z 

z 

2 

SO(3) 

0 

6 

Z 2 

0 

z 

3 

SO(4) 

U(2) 

6 

Z 2 

z 

z 

2 

SO(5) 

SO(3) X SO(2) X Z 2 

14 

Z 2 

z 

Z 2 

6 


Table 1: Models. We list the models considered in this paper, in which a global symmetry 
gronp G is spontaneonsly broken to a snbgronp H. N is the nnmber of real scalar helds; 
the twelve real helds in the SU(3) theory are arranged into two complex three-vectors, 
while the other representations are explicitly real. The lower homotopy gronps of the 
vacnnm manifold G/H are listed, as well as the dimensionality of G/H (corresponding 
to the nnmber of massless Goldstone bosons). The models are gronped into two classes, 
depending on whether the model snpports strings [ki{G/H)] and/or monopoles [7r2{G/H)]. 
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Figure Captions 


Figure One. These plots represent energy density contonrs in the x-y plane of oblate 
conhgnrations in the SO(4)/SO(3) and SO(3)/SO(2) theories. The textnres are obtained 
from those with spherically symmetric energy densities by shrinking the conhgnrations 
by 50% along the x-axis. The energy density in the SO(4)/SO(3) is simply deformed 
along with the coordinates, while in the SO(3)/SO(2) model the maxim um energy density 
describes a ring in the y-z plane. In these hgnres, the three-dimensional conhgnration is 
obtained by rotation aronnd the x-axis. 

Figure Two. The integrated energy densities dnring the evolntion of textnres with ini¬ 
tially spherically symmetric energy densities in the fonr theories which do not snpport 
strings or monopoles. Each form of energy density is integrated over the volnme of the 
simnlation region; the total energy is the snm of the integrated potential, kinetic and 
gradient energies. The horizontal axis is measnred in timesteps of the simnlations. 

Figure Three. The integrated energy densities dnring the evolntion of textnres with 
initially oblate conhgnrations (obtained from the symmetric conhgnrations by deforming 
along one axis) in the fonr theories which do not snpport strings or monopoles. 

Figure Four. Qnadrnpole moments of the total energy density for collapse of initially 
oblate conhgnrations in the fonr theories which do not snpport strings or monopoles. 

Figure Five. The integrated energy densities dnring the evolntion of textnres with ini¬ 
tially spherically symmetric energy densities in the fonr theories which do snpport either 
strings or monopoles. (These textnres are referred to as “symmetric”, althongh only in 
SO(5)/[SO(3) X SO(2) X Z 2 ] is the conhgnration trnly spherically symmetric.) 

Figure Six. Qnadrnpole moments of the total energy density for collapse of conhgnrations 
with initially symmetric energy densities in the fonr theories which do snpport either strings 
or monopoles. Note that nonzero qnadrnpoles develop for three of the fonr models. 
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Figure Seven. Three-dimensional contours of total energy density and potential energy 
at one moment during the collapse of the SO(3)/SO(2) conhguration with an initially 
symmetric energy density. The shaded regions are those with large potential energies, 
representing monopoles and anti-monopoles. The lower total-energy contour describes a 
prolate spheroid, while the higher contour describes two roughly hemispherical regions 
associated with the inner mo nopole/anti-mo nopole pair. 

Figure Eight. Three-dimensional contours of total energy density and potential energy at 
one moment during the collapse of the SO(3)/0 conhguration with an initially symmetric 
energy density. The shaded region is one of large potential energy, representing a loop of 
cosmic string. The total-energy contours describe roughly spheroidal regions. 

Figure Nine. The integrated energy densities during the evolution of textures with 
initially oblate energy densities (obtained from the symmetric conhgurations by deforming 
along one axis) in the four theories which do support either strings or monopoles. 

Figure Ten. Quadrupole moments of the total energy density for collapse of initially 
oblate conhgurations in the four theories which do support either strings or monopoles. 
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